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CNJ | Abstract 

We investigate the effects of very strong magnetic fields upon the equation of state 
\q • of dense bayonic matter in which hyperons are present. In the presence of a magnetic 

field, the equation of state above nuclear density is significantly affected both by 
Landau quantization and magnetic moment interactions, but only for field strengths 
B > 5 x 10 18 G. The former tends to soften the EOS and increase proton and 
lepton abundances, while the latter produces an overall stiffening of the EOS. Each 
results in a supression of hyperons relative to the field-free case. The structure of 
a neutron star is, however, primarily determined by the magnetic field stress. We 
utilize existing general relativistic magneto-hydrostatic calculations to demonstrate 
that maximum average fields within a stable neutron are limited to values B < 
1— 3x 10 18 G. This is not large enough to significantly influence particle compositions 
or the matter pressure, unless fluctuations dominate the average field strengths in 
the interior or configurations with significantly larger field gradients are considered. 
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Recent discoveries that link soft 7-ray repeaters and perhaps some anomalous 
X-ray pulsars with neutron stars having ultrastrong magnetic fields - the so 
called magnetars (see Table 1 in [1] for a summary) - have spurred theoretical 
studies of the effects ultrastrong magnetic fields have on the equation of state 
(EOS) of neutron-star matter and on the structure of neutron stars (cf. [1-3] 
and references therein) . The magnetic field strength B needed to dramatically 
affect neutron star structure directly can be estimated with a dimensional 
analysis [4] equating the magnetic field energy E B ~ (47i\R 3 /3)(.E> 2 /87r) with 
the gravitational binding energy E B E _ ~ 3GM 2 /5R, yielding 

B^xW^-V-^fV (1) 
V1.4 M ; UO km/ ' v ; 

where M and R are, respectively, the neutron star mass and radius. 

To date, studies of the effects of such ultrastrong magnetic fields on the EOS 
and on the structure of neutron stars have been largely limited to cases in 
which the strongly interacting component of matter is comprised of nucle- 
ons [3,5-7]. Our objective in this paper is to investigate the effects of magnetic 
fields on the EOS of matter containing strangeness-bearing hyperons. Towards 
this end, we utilize the theoretical formalism for the EOS developed in Ref. [3]. 
This allows us to consistently incorporate the effects of magnetic fields on the 
EOS in multicomponent, interacting matter. For the first time, the anoma- 
lous magnetic moments of baryons is included, using a covariant description. 
A further objective is to explore the consequences of strong magnetic fields 
for neutron star structure and maximum masses, which is accomplished using 
the results of general relativistic structural calulations in Ref. [1]. 

Charge neutral, beta-equilibrated, neutron-star matter contains both nega- 
tively charged leptons (e and ji) and positively charged baryons (p and, at 
higher densities, S + ). Magnetic fields quantize the orbital motion (Landau 
quantization) of these charged particles. When the Fermi energy of the p or 
S + becomes significantly altered by the magnetic field, the pressure and com- 
position of matter in beta equilibrum are affected. Landau quantzation has a 
significant effect when B* = B / B e c ~ 10 5 (B e c = hc/(e A e 2 ) = 4.414 x 10 13 G 
is the critical electron field) [3,5-7]. Higher fields lead to a general softening 
of the EOS relative to the case in which magnetic fields are absent. 

In strong magnetic fields, contributions from the anomalous magnetic mo- 
ments of the baryons (denoted hereafter by n b ) must also be considered. 
The measured magnetic moments ji^ for nucleons and strangeness-bearing 
hyperons are given in Table 1. With few exceptions, the anomalous mag- 
netic moments are similar in magnitude. As discussed below, the energies 
\n n B\ ~ 2.7 x \Q- A B* MeV and \k p B\ ~ 2.5 x 10~ 4 B* MeV measure the 
changes to the nucleon Fermi energies, and \n n + k p \B ~ 1.67 x 1CT 5 £?* MeV 
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measures the change to the beta equilibrium condition. Since the Fermi ener- 
gies exceed 20 MeV for the densities of interest, it is clear that contributions 
from the anomalous magnetic moments become significant for B* > 10 5 . For 
such large fields, complete spin polarization of the neutrons occurs, producing 
an increase in degeneracy pressure and an overall stiffening of the EOS that 
overwhelms the softening induced by Landau quantization [3]. 

Similary, in the presence of strong magnetic fields, all of the hyperons are sus- 
ceptible to spin polarization due to magnetic moment interactions with the 
field. The net result of the opposing effects of degeneracy and Landau quan- 
tization will depend on the relative concentrations of the various hyperons, 
which in turn depend sensitively on the hyperon-meson interactions for which 
only a modest amount of guidance is available [9-11]. It is one of the purposes 
of this work to ascertain the influence of feedback effects due to mass and 
energy shifts in multi-component matter. 



TABLE 1: The measured magnetic moments of spin— ^ baryons from Ref. [8]. 
fiN = eh/(2m p ) = 3.15 x 10 -18 MeV G _1 is the nuclear magneton and the 
anomalous magnetic moment Kb = (f^b/f^N — Qb m p/ m b) ^N- 



Species 
b 


Mass 
(MeV) 


Charge 
Qb 


Magnetic Moment 


Anomalous Moment 


P 


938.3 


1 


2.79 


1.79 


n 


939.6 





-1.91 


-1.91 


A 


1115.7 





-0.61 


-0.61 


£+ 


1189.4 


1 


2.46 


1.67 


S° 


1192.6 





1.61 


1.61 


s~ 


1197.4 


-1 


-1.16 


-0.38 




1314.9 





-1.25 


-1.25 




1321.3 


-1 


-0.65 


0.06 



To describe the EOS of neutron-star matter, we employ a field theoretical ap- 
proach (at the mean field level) in which the baryons interact via the exchange 
of a-u-p mesonsQ We consider two classes of models that differ in their high 
density behavior. The two Lagrangian densities are given by [12,13] 



The qualitative effects of strong magnetic fields found in the field-theoretical 
models also exist in non-relativistic potential models. This is because the phase 
space of charged particles is similarly affected in both approaches by magnetic fields. 
The effects due to the anomalous magnetic moments would, however, enter linearly 
in a non-relativistic approach, and would thus be more dramatic in that case. 
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V m b J 

Cii = ( 1 + — ) C b - %m b V b + C m + C t ■ (2) 

V 7716/ 

The baryon, lepton, and meson Lagrangians are given by 

C b = ^ b (i^d» + q bltl A» - g^u" ~ 9pbT 3b l^ - W^F^) f b , 
C e = ^ (i^ff* + qa^A") fa , and 

C m = \d»ad»a - l -mla 2 - U(a) + \m 2 ^^ - 

+ \rn 2 p p,ff - - \F^F, U , (3) 

where ty b and ipe are the baryon and lepton Dirac fields, respectively. The 
index b runs over the baryons n, p, A, E _ , S°, £ + , S _ , and S°. (Neglect of 
the fl~ and the A quartet, which appear only at very high densities, does not 
qualitatively affect our conclusions.) The index £ runs over the electron and 
muon. The baryon mass and the isospin projection are denoted by m b and r 3b , 
respectively. The mesonic and electromagnetic field tensors are given by their 
usual expressions: = d fl u u - d u u^, = d^p u - d u p^ and F^ = d^A u - 
dyA^. The strong interaction couplings are denoted by g, the electromagnetic 
couplings by q, and the meson masses by m all with appropriate subscripts. 
The quantity U(a) denotes scalar self-interactions and is taken to be of the 
form U(a) = (bm n /3)(g aN a) 3 + (c/±)(g aN a) A [12,14,15]. 

The anomalous magnetic moments are introduced via the minimal coupling of 
the baryons to the electromagnetic field tensor for each baryon. Their contri- 
butions in the Lagrangian are contained in the term — K b ^ b a ^F^^ b , where 
Qiiv = I hV 4 ' ^"l - Although the electromagnetic field is included in Ci and Cn, 
it assumed to be externally generated and only frozen-field configurations will 
be considered. 

The meson- nucleon <7ijv(z = cr, oj, p), and scalar self-interaction couplings, b and 
c, were chosen to reproduce the binding energy B/A, the nuclear saturation 
density n s , the Dirac effective mass M*, and the symmetry energy a s . The 
values of the couplings used in the models considered are shown in Table 2, 
and the associated nuclear matter properties used to calculate these couplings 
are shown in Table 3. 

The meson-hyperon couplings are assumed to be fixed fractions of the meson- 
nucleon couplings, i.e., gm = x^giN, where for each meson %, the values of 
XiH are assumed equal for all hyperons H . The values of xm are chosen to 
reproduce the binding energy of the A at nuclear saturation as suggested by 
Glendenning and Moszkowski (GM) [9] and are also given in Table 2. The 
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TABLE 2: Nucleon-meson and hyperon-meson coupling constants for the GM1-3 
and ZM models. The hyperon to nucleon coupling ratios were taken from Ref. [16]. 



Model 


9aN/m a 
(fm) 


9u>N/mu> 
(fm) 


g P N/in p 
(fm) 






XpH 


b 


c 


GM1 


3.434 


2.674 


2.100 


0.6 


0.653 


0.6 


0.002947 


-0.001070 


GM2 


3.025 


2.195 


2.189 


0.6 


0.659 


0.6 


0.003478 


0.01328 


GM3 


3.151 


2.195 


2.189 


0.6 


0.659 


0.6 


0.008659 


-0.002421 


ZM 


2.736 


1.617 


2.185 


1.0 


1.0 


1.0 


0.0 


0.0 



TABLE 3: Nuclear matter properties used to fit coupling constants for the GM1-3 
Glendenning and Moszkowski (GM1-3) Zimanyi and Moszkowski (ZM) models. 



Model 


n s 
(fm" 3 ) 


-B/A 
(MeV) 


M*/M 


K 
(MeV) 


(MeV) 


GM1 


0.153 


16.30 


0.70 


300 


32.5 


GM2 


0.153 


16.30 


0.78 


300 


32.5 


GM3 


0.153 


16.30 


0.78 


240 


32.5 


ZM 


0.160 


16.00 


0.86 


225 


32.5 



high-density behavior of the EOS is sensitive to the strength of the meson 
couplings employed [10] and the models chosen encompass a fairly wide range 
of variation. Models GM1-3 employ linear scalar couplings (£/), while the 
Zimanyi and Moszkowski (ZM) model employs a nonlinear scalar coupling 
(£//), which is reflected in the high density behaviors of m*/m n . Comparison 
of results from the GM1-3 and ZM models allows us to contrast the effects of 
the underlying EOS. 

The field equations (not explicitly displayed here) follow from a natural ex- 
tension of those used in Ref. [3] to the case in which hyperons are present. In 
charge neutral and beta equilibrated matter, the conditions 

1bn b + qm = and \i h = B b fi n - q h \i e , (4) 
b i 

where is the baryon chemical potential (not to be confused with the mag- 
netic moment //& of Table 1) and is the baryon number, also apply. The 
energy spectra, number densities, scalar number densities, energy densities, 
and total pressure are given by straightforward generalizations of their ana- 
logues in Ref. [3]. 

In Fig. 1, results are shown for physical quantities of interest for the baseline 
model GM3, when the effects of Landau quantization are included but those of 
anomalous magnetic moments are ignored. As the magnetic field is varied, the 
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u = n b /n u = n b /n 

Fig. 1. Matter pressure P m , neutron Dirac effective mass m*/m n , and concentrations 
Yi as functions of density (no = 0.16 fm~ 3 is the fiducial nuclear saturation density) 
for the model GM3. Each is shown for magnetic field strengths of B* = 0, 10 5 , and 
10 6 . Pf = B 2 /8tt is the pressure due to magnetic field stress for the value B* = 10 5 . 
The inset shows P m as a function of the matter energy density e m . In the right 
panels results for leptons are shown by dashed lines while those for baryons are 
shown by solid lines. 

generic behavior of the matter pressure P m , effective neutron mass to*, and 
particle abundances are the same as those found in Ref. [3], in which matter 
containing only nucleons was considered. When hyperons appear, the pressure 
becomes smaller than the case of pure nucleons for all fields. The effective 
nucleon mass decreases slightly when hyperons appear, again independently 
of the field strength. The appearance of hyperons produces pronounced effects 
on the nucleon and lepton abundances. 
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The cases B* = and B* = 10 5 are nearly indistinguishable, except for 
densities less than nuclear density (u < 1), where "cavitation" due to Landau 
quantization occurs. This effect has been discussed in Ref. [3] and earlier work 
(cf., Refs. [5,6]). However, the structure of a magnetized neutron star will be 
mostly affected by contributions from the magnetic field stress, Pf = B 2 /8n = 
4.814 x 10 _8 -B* MeV fm~ 3 , which greatly exceeds the matter pressure P m at 
all relevant densities for B* > 10 5 , as shown in Fig. 1 for B* = 10 5 . 

Hyperons appear near twice nuclear density for B* < 10 5 in this model. For 
B* > 10 5 (B > 5 x 10 18 G), the same value needed to noticeably affect the 
matter pressure and nucleon abundances, significant effects on hyperon thresh- 
old densities and abundances are apparent. Landau quantization increases the 
proton abundance, and hence the electron concentration due to charge neu- 
trality. Therefore /i e — fj, n — ji p is increased. However, at fixed baryon density, 
the neutron abundance must then decrease, which significantly lowers ji n to 
a greater extent than fi n — fi p is increased. The net effect is a suppression of 
hyperons. In particular, note that // s - = 2/i n — fi p and //a = Un- 
it should also be noted that, compared to the case without hyperons (cf. 
Ref. [3]), the pressure is reduced for all field strengths, as is the nucleon effec- 
tive mass. The magnitude of the nucleon effective mass can be important in 
determining the chemical potentials, as shown below. 

We have verified that the qualitative effects of large fields on the composition, 
presure and effective mass are similar for the other three models considered, in 
the case that only Landau quantization is considered. The relative suppression 
of individual hyperons varies, however; for example in models GM2 and ZM 
the A and £~ threshold densities reverse at large fields relative to GM1 and 
GM3. Importantly, however, the value B* ~ 10 5 represents a threshold for 
field effects to become significant, independently of the EOS. 

The inclusion of anomalous magnetic moments produces even larger suppres- 
sions of hyperons, as displayed in Fig. 2 for the baseline model GM3. One 
can qualitatively understand this trend by considering the effects of magnetic 
moments on the baryon chemical potentials, and then examining the effects on 
the threshold of the first hyperon to appear. At low temperatures, the chem- 
ical potentials can be identified with the Fermi levels of the baryonic energy 
spectra. The energy spectrum for charged baryons is given by [3] 



where the quantity z/ = n + s/2 + l/2 characterizes the so-called Landau level 




= y k\ + ^ml 2 + 2uq b B + sk p B) + g Wp u° + t 3b g Pp p° , 
k 2 z vq b B n , n 

- m b + tt^ + — ~ + SK b B + + h h 9 Pb p , 



(5) 
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02468 10 02468 10 



u = n b /n u = n b /n 

Fig. 2. Same as Fig. 1, but with magnetic moment interactions included. 

with n and s (±1) being the principle quantum number and "spin" quantum 
number respectively, and t 3b is the isospin projection of baryon b. The quanti- 
ties u; and p° are mean field solutions of the field equations corresponding to 
Eq. (3) (cf. Ref. [3]). The rightmost expression is valid in the nonrelativistic 
limits kl <^ml , 2v\qi\B <C ml and (k^B) 2 <C ml , which are appropriate for 
the conditions of interest. Similarly the energy spectrum of neutral baryons 
is given by 



E b ,s = sjk 2 z + (yfmi 2 + k 2 x + k 2 + sK b B) + g„ b u° + t 3b g Pb p° , 
k 2 

~ m b + ^* + SKbB + 9 ^ u ° + tMnP • ( 6 ) 
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For both charged and neutral baryons, the chemical potentials of nearly de- 
generate systems are set by the lowest energy spin state. The contribution to 
neutral baryon energies, and to their chemical potentials, due to the anoma- 
lous magnetic moment is therefore — \ki,\B in Eq. (6). The situation of charged 
baryons is more complicated, because of the additional contribution due to the 
Dirac moment. However, given that the value of n from Landau quantization 
generally exceeds |/c&|///&, the shift in charged baryon chemical potentials is 
not as significant as it is for neutral baryons. In matter containing only nucle- 
ons, the predominant effect of the anomalous magnetic moments is therefore 
to decrease /i n , and thus the neutron abundance, relative to the case in which 
the anomalous moments are ignored (compare Figs. 1 and 2). A further con- 
sequence is that the threshold density for the first hyperon, S~, is increased 
since //£- = 2/i n — fi p . Although the chemical potential of the S _ also receives 
a contribution from its anomalous magnetic moment, the magnitude of the 
change is small compared to that due to the neutron. Note again that field 
strengths B* > 10 5 are necessary to produce appreciable effects. We have ver- 
ified that these conclusions hold for the other three EOSs we considered for 
the case when magnetic moments are included. 

While most hyperon threshold densities are increased by the inclusion of 
anomalous magnetic moments, those of S° and S + , the latter to a lesser 
extent, are decreased for B* > 10 5 . The complex interplay that exists when 
several baryons are simultaneously present makes a simple explanation diffi- 
cult, however. 

Finally, as noted in Ref. [3], the inclusion of magnetic moments produces 
an increase in pressure and effective mass at large fields (B* > 10 5 ), due 
to increased baryon degeneracies [3]. These increases more than offset the 
reductions induced by the inclusion of hyperons. 

Still to be addressed, however, is whether or not stable stellar configurations 
can exist in which the magnetic field is large enough (B > 5 x 10 18 G) that 
the properties of matter are significantly affected by the magnetic field. One 
step in this direction was recently undertaken in Ref. [1], in which the lim- 
its of hydrostatic equilibrium for axially-symmetric magnetic fields in general 
relativistic configurations were analyzed. As discussed in detail in Ref. [2], in 
axially symmetric field configurations with a constant current function, the 
magnetic field contributes a centrifugal-like contribution to the total stress 
tensor. This can be understood by noting that for the artificial field geometry 
considered (namely B oc z), a superconducting fluid can move along field lines 
but not across them. Thus the "pressure" associated with the magnetic field 
will only act equatorially and not vertically. This flattens an otherwise spher- 
ical star, and for large enough fields, decreases the central (energy) density 
even as the mass is increased. For large enough fields, the star's stability is 
eventually compromised. 



9 




Fig. 3. Limits to hydrostatic configurations for neutron stars permeated by axially 
symmetric magnetic fields. The upper (lower) panel is for the EOS GM3 including 
(excluding) hyperons. In each panel, the lower heavy solid curve is the standard 
mass-radius relation for field- free stars. The upper heavy solid curve represents the 
largest gravitational mass possible for a given equatorial radius as the magnetic 
fields are increased, for the indicated current function /. The heavy dashed curve 
is the locus of minima of the contours of fixed baryon mass configurations: Ref. [1] 
suggested that this could be the limit to dynamical stability. Thin solid lines are 
contours of the maximum magnetic field strength in units of 10 18 G. Similarly, thin 
dashed lines are contours of maximum mass-energy density in units of the energy 
density po at the nuclear saturation density. 

Hydrostatically stable configurations (of which some may not be stable to 
dynamical perturbations) for the GM3 EOS, both excluding and including 
hyperons, are shown in mass-equatorial radius space in Fig. 3. The lower 
thick black line in each panel is the usual field-free, spherical result for the 
mass-radius relation. The upper thick black line represents the largest possi- 
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ble stable mass for a given equatorial radius as the internal field is increased. 
Large axially-symmetric fields tend to yield flattened configurations, and if 
large enough, shift the maximum densities off-center in the shaded regions in 
the figure. This results in toroidal shapes with low-density centers. As dis- 
cussed in Ref. [1], regions to the left of the thick dashed line in this figure 
are likely unstable to small amplitude perturbations and hence unrealistic as 
stable physical configurations. 

Superimposed on this figure are contours of maximum mass-energy density 
and magnetic field strength. For the GM3 model, hyperons appear at zero 
field at about twice nuclear density. Therefore, the portions of the two panels 
in which the maximum density is below about 3 times nuclear density are 
nearly identical. For the particular choice of a constant current function, but 
independently of the EOS, it is found that the maximum value of the magnetic 
field within the star cannot exceed about 3 x 10 18 G, which corresponds to 
5*~7x 10 4 . In the case of GM3, the maximum value is B ~ 1.8 x 10 18 G. 
This field strength is not nearly large enough to produce appreciable effects 
on hyperon or nucleon compositions, with or without the inclusion of anoma- 
lous magnetic moments. The dominant effect of the field arises through the 
magnetic field stress, which effectively dominates the matter pressure below 
a few times nuclear saturation density, depending on the field's orientation. 
Whether other choices of current functions, or the relaxation of the condition 
of axial symmetry, will alter these conclusions is not clear, as the answer will 
likely depend on the detailed nature of the field configurations considered. 

To date only a few cases of spatially varying current functions have been 
explored [2], but these have all been limited to axial symmetry. Furthermore, 
the cases explored have relatively small spatial gradients so that the ratio of the 
maximum field to the average field within the star is not large. It might be that 
the shapes of the stars will significantly change with a different field geometry. 
It is even possible to imagine a disordered field where (B 2 ) is significantly 
larger than (B) 2 . In this case the pressure will be dominated by fluctuations 
in the field, but the stars will tend to be spherical. It is possible that strong 
magnetic fields may be held in the core for periods much longer than the ohmic 
diffusion time due to interactions between the magnetic flux tubes and the 
vortex tubes expected to be present in a superconducting, superfluid, rotating 
neutron star [17]. Therefore, although the results to date imply that average 
fields within a neutron star cannot exceed 1 — 3 x 10 18 G before stability is 
compromised, the value of the maximum possible field at any point within the 
star is undetermined. As a result, it is possible that the effects of magnetic 
fields upon the equation of state or particle composition will be important in 
some magnetic field geometries. 
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87ER40317 (for JML) are gratefully acknowledged. We are grateful to Chris- 
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